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I.   INTRODUCTION 

During  the  last  fifteen  years,  there  have  been  many 
advancements  in  the  development  of  highly  effective  tacti- 
cal weapons  for  use  at  intermediate  range  and  with  high 
average  flight  speed.   One  of  the  simplest  engines/motors 
for  such  use  is  the  Solid  Fuel  Ram  Jet  (SFRJ) . 

The  SFRJ  most  often  consists  of  a  solid  fuel  grain  which 
provides  the  walls  for  the  combustion  grain  [Ref .  1] .   A 
sudden  expansion  at  the  air  inlet  of  the  combustor  can  be 
used  to  provide  flame  stabilization  by  inducing  recircu- 
lation into  the  air  flow  (Fig.  1.11.   Combustion  can  be  sus- 
tained throughout  the  grain  provided  that  the  step  height 
(h)  is  sufficiently  large.   However,  the  larger  the  step 
height,  the  larger  the  losses  of  the  inlet  stagnation  pressure 

Due  to  the  sudden  expansion,  at  least  two  distinct  flow 
fields  are  generated  within  the  fuel  grain  (Fig.  1.2).   In 
the  first  field,  the  recirculation  zone,  the  flow  is  highly 
turbulent  and  usually  fuel  rich.   This  hot  gas  region  pro- 
vides the  energy  necessary  to  sustain  the  combustion  proc- 
ess which  must  occur  further  downstream.   Downstream  of  the 
flow  reattachment,  a  turbulent  boundary  layer  develops  and 
the  combustion  is  diffusion  controlled.   A  diffusion  flame 
emanates  from  the  recirculation  zone  and  remains  within 
the  developing  turbulent  boundary  layer  between  the 
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Figure  1.2.   Effects  of  Sudden  Enlargements  on  a  SFRJ 

fuel-rich  zone  near  the  wall  and  the  oxygen-rich  zone  near 
the  central  core. 

There  have  been  continuing  research  efforts  concerning 
the  behavior  of  the  SFRJ  at  the  Naval  Postgraduate  School. 
Both  mathematical  and  experimental  efforts  have  taken  place 
to  determine  the  effects  of  design  as  well  as  operational 
variables  of  the  system  to  be  obtained. 
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II.   DESCRIPTION  OF  THE  COMPUTER  MODEL 


A.  INTRODUCTION 

At  the  NPS  there  have  been  two  basic  computer  codes 
used  to  model  the  flow  within  solid  fuel  ramjets.   The 
first  one  was  based  on  stream  function  (.ty)    and  vorticity  (co) 
[Ref.  3] .   It  predicted  with  reasonable  accuracy  the  velocity 
and  temperature  distributions  but  did  not  predict  pressure 
distributions  accurately  and  was  not  easily  extended  to 
complex  geometries  or  to  three-dimensional  flows.   The  second 
computer  model  is  based  on  the  "primitive"  variables,  pres- 
sure (p)  and  velocities  (u,v)  [Ref.  2J .   It  has  also  been 
used  to  predict  the  effects  of  fuel  properties  on  the  SFRJ 
performance  and  to  evaluate  the  effects  of  different 
geometries  as  well  as  operating  conditions. 

B.  OVERVIEW  OF  THE  PRIMITIVE  VARIABLE  PROGRAM 
1 .   Assumptions 

The  CHAMP  ION/ 2/E/F IX   computer  program  developed  by 
Pun  and  Spalding  [Ref.  2]    has  been  used  as  the  basis  for  the 
primitive  variable  model.   The  flow  has  been  assumed  to  be 
steady  (solid  fuel  grain  dissipation  is  negligible! ,    sub- 
sonic, recirculating,  axisymmetric  and  to  have  constant 
specific  heats.   In  order  to  calculate  the  effective  viscosity, 
a  modified  Jones-Launder  [Ref.  4 J  two-parameter  turbulence 
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model  has  been  used.   It  uses  five  empirical  constants 
(Table  II)  and  requires  that  two  additional  variables 
(turbulent  kinetic  energy  (k)  and  turbulence  dissipation 
rate  (e))  be  evaluated.   The  effective  viscosity  was  calcu- 
lated using  the  formulas 


V  Jin        ~        Ui     +  Uj.  (2.1) 

Heff     Mlam   Kt.. 


where 


ufc   =  CD  p  k2/z  C2.2) 


TABLE  I 
k-e  Turbulence  Model  Empirical  Constants 


c,       c        c_       a,  j.j.  o        ~j- 

1        2        D        k,eff        e,eff 


1.43      1.92       .09        1.0  1.3 


The  combustion  was  considered  to  be  mixing- limited 
(infinitely  fast  kinetics)  with  a  simple  one-step  chemical 
reaction  of  the  form: 

1  Kg  fuel  +  i  Kg  oxidizer  •*■      (l+i)kg  products     (2.3). 

Four  species  were  concerned:   oxygen,  nitrogen,  fuel 
and  products.   Because  of  the  infinitely  fast  kinetics 
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Governing  Equation  Parameters 
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assumed,  -  no  oxygen  could  exist  at  the  fuel  surface  and  at 
the  surface  was  assumed  to  be  isothermal.   The  turbulent 
Prandtl  and  Schmidt  numbers  were  taken  to  be  equal  to  unity 
and  therefore  the  turbulent  Lewis  number  and  laminar  Prandtl 
number  were  also  taken  to  be  equal  to  unity. 

C.   GOVERNING  EQUATIONS 

With  the  above  assumptions,  the  governing  equations  for 
axi symmetrical  flows  can  be  written  as  follows: 


f-(pucj>)  +  i  f-(prv4>)  -  f-(I\  !$)-£■  f-(rr.  |£)  +  S  .  =  0   (2.4) 
dx        r  dr         dx   9  dx    r  3r    9  dr     9 


where  9  stands  for  the  independent  variables  u,  v,  k,  e, 

h,  itu  -m  /i  and  hlt  .   The  first  two  terms  are  the  "convec- 
fu   ox'       TSI2 

tion  terms",  the  next  two  the  "diffusion  terms".   r,  is  the 
effective  exchange  coefficient  for  turbulent  flow  (Table  I) 
and  S   is  the  "source  term"  (Table  II).  .   The  stagnation 
enthalpy  for  unity  Prandtl  and  Schmidt  numbers  is 


h  =   h  +  (u2+v2)/2  +  K  C2.5). 


where  for  nonreacting  flows 


h   E   c  T  C2.6) 

P 


and  for  reacting  flows 
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AH 

hEm   =£■  +  c  (T  -  T   -)  (2.7) 

ox  1     p      ref  v    ' 


Temperature  has  been  calculated  using  Eqs .  (2.5) 
and  (2.7).   Density  has  been  calculated  using  the  perfect 
gas  law 


=  &  (2-8) 


Since  the  quantities  h  (stagnation  enthalpy) ,  m^ 

N2 

and  x  =  m«   -  m  /i 

A    fu    ox' 

have  identical  governing  differential  equations  and,  in 
appropriate  dimensionless  form,  identical  boundary  condi- 
tions, a  representative  equation  had  to  be  solved  only  once 
The  dimensionless  enthalpy  was  selected  to  be 


h.   -  h 
H  =   ^±2 £2.9) 

h.   -  h - 

in    fg 


where  h.      is  the  stagnation  enthalpy  deep  in  the  fuel  grain 
and  h.   is  the  stagnation  enthalpy  of  the  air  at  the  inlet 
dump  plane.   The  corresponding  dimensionless  forms  for  the 
other  quantities  were: 


-    u     __2in 2_  i2 

2m      2fg 

and 
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x      = 

x    -    X. 

in 

x,-      -   x . 

fg           in 

fu 

= 

X 

m 
ox 

=      0 

f  u 

= 

0 

m 
ox 

=      -x. 

1 

(2.11) 


for   x  >  0 

(2.12) 
for   x  <  0 


1 .   Boundary  Conditions  and  Solution  Procedure 

Fixed  boundary  conditions  (inlet  velocity  profile, 
et  al . )  were  specified  at  the  desired  or  experimentally 
determined  values.   Specified  gradient  boundary  conditions 
were  handled  by  setting  the  appropriate  convection/diffusion 
coefficients  to  zero  in  the  finite  difference  equation 
("breaking  the  link")  and  then  entering  the  appropriate 
gradient  through  linearized  "false"  source  terms.   The. 
geometry  as  well  as  the  appropriate  boundary  conditions  are 
summarized  in  Figure  (2.1). 

"Plug  flow"  was  assumed  at  the  air  inlet,  although 
this  was  not  a  computer  program  limitation.   The  turbulence 
kinetic  energy  was  assumed  to  be  uniform,  with  a  value  which 
corresponded  to  the  approximate  turbulence  intensity  of  the 
inlet  flow.   Radial  and  axial  gradients  were  set  equal  to 
zero  on  the  centerline  and  at  the  exit,  respectively.   All 
the  nonreacting  solid  boundaries  were  considered  adiabatic 
with  both  u  and  v  velocity  components  equal  to  zero  (no-slip 
condition) . 
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For  simplicity,  the  boundary  layer  was  assumed  to  be 
composed  of  a  laminar  sublayer  which  is  then  developed  to 
turbulent  without  assuming  a  "buffer  zone".   For  the 
laminar  sublayer, 


u+  =  y+  (2.13) 


P     "P 


where 


u  /t   /py 

+                    P  ^                 +_w/Kjrp 

u        =      *-—  and  y        =      *- 

p    /T7p  p 


and  for  the  rest  of  the  layer 

£n (Ey ' ) 


+  , 


u   =  2-  (2.14) 

P  K 


For  E  =  9.0  and  k  =  0.4  (values  found  from  experimental 
data) 


u+   =   5.5  +  2.5  in   y+  (2.151 

P  P 


Solving  the  above  two  equations,  it  is  found  that 

transition  from  laminar  to  turbulent  flow  occurs  at  y  =11.5 

The  value  of  y  at  each  near-wall  node  is  found  from  the 

^P 

relation: 


5   t   1/2 
y+   =   -£—{-£)  (2.16) 

P     ^lam  p 
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For  y   <  11.5 
P 


y,   -u 
=   _lam^ 

wo 


and  for  y   >  11.5 
P 


t    =  ci/2  p  k  (2.18) 

w       D    K   p 


The  value  of  t   was  assumed  to  be  uniform  from  the  wall  to 

w 


the  near-wall  grid  point.   So, 


+      1/4     1/2   6  ,0  1QN 

Yp  =  C°   P  Kp   ^  l2-19) 


Then,  the  shear  stress  is  calculated  using  the  formula 
[Ref.  3] 


u 

w       D    ""  "p         D     p 


-   cJ/2  p  <    =   p  =V4  Kl/2(^ 


u 

p 

.  _     1/4  -1/2 
,  /9     E  p  c  '   k  ' 

=   Kpu   k  '   £n( 2 — )      (2.20) 

P   P  ylam 


where,  as  noted  previously, 


£n  (Ey+) 
u   =   1- 


P  K 

Due  to  the  steep  gradients  of  the  properties  in 
turbulent  flows  near  the  solid  boundaries,  the  source  terms 
for  k  and  e  at  the  near-wall  nodes  were  expressed  in  terms 
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of  the  wall  shear  stress.   Also,  t   provides  the  boundary 

condition  for  the  u  and  v  equations. 

The  turbulence  dissipation  rate  e   at  a  near-wall 

P 

node  p  is  given  by  the  relation  [Ref .  3] 


E    =   c3/4  k3/4/k5   =   k3/2/2.436  (2.21) 

p       D     p  p 


where  the  length  scale  is  presumed  to  be  proportional  to  the 
incremental  distance  from  the  wall, 


I      =   c^/4  k  6  (2.22) 


It  was  found  that  when  using  a  sudden  expansion 

geometry  in  reacting  flows  the  near-wall  dissipation  rate 

3/2 
had  to  be  increased  on  the  step  face  (e   =  k  '    /0.46)  and  the 

P    P 
grid  spacing  adjacent  to  the  grain  had  to  be  fine  (y   <  11.5) 

in  order  to  obtain  a  temperature  distribution  in  qualitative 
and  quantitative  agreement  with  experimental  data  [Ref.  1J . 
Equation  C2.17)  implies  that  the  wall  shear  stress  is  calcu- 
lated assuming  a  linear  velocity  profile  when  y   <  11.5. 
Therefore,  a  near-wall  grid  point  can  lie  within  the  laminar 
sublayer  but  the  source  terms  for  k  and  e  imply  that  y  f-p/u-i  m 
is  much  greater  than  unity.   This  precludes  y   from  being 
significantly  less  than  11.5. 

For  reacting  flows,  the  boundary  conditions  for  the 
dimensionless  properties  (Eqs.  (2.9)  through  (2.12))  were  zero 
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at  the  inlet  and  unity  "deep"  in  the  fuel  grain.   These 
properties  were  considered  to  have  zero  gradients  on  non- 
reacting  surfaces. 

The  assumptions  made  for  reacting  flows  (unity 
Prandtl  and  Schmidt  numbers,  simple  chemical  reaction,  con- 
stant specific  heat  and  stagnation  enthalpy  defined  in  Eq . 
(2.5))  result  in  a  general  boundary  condition  for  all 
"conserved"  properties  (A  )  on  a  surface  which  has  mass 
transfer  given  by  the  relation  [Ref .  2] 

mT         =   (T,  x-=-)u  /(A   u   -  <j>  c    )  (2.23) 

Tdw       A  8r   bw7  Yc,bw   Yc,fg 

where  <p     may  represent  h/itLj   or  x  =  mf   -  m  /i. 

A  mass  transfer  conductance  (g)  is  often  defined 
from  the.  relation 


3<t>c 

(r.  t^)k   =  g(<i>   -  <(i  .J  (2.24) 

<f>  9r   bw     3  Yc,°°    Tc,bw 


where  $   is  defined  as  the  free-stream  value  of  the  property 

Substituting  Eq.  (2.24)  in  Eq.  (2.23),  the  mass 
transfer  from  the  blowing  wall  is 

m£    =   g(*    "  <f>  u  )/(*   u   -  <J>   *  )  (2.25) 

bw     a  Tcp    Ycbw"  Yc,bw   Yc,fg 

where 
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BP  =  *cp  -  ♦c.lw  (226) 

^Cjbw    c,fg 


represents  the  mass  transfer  parameter  (or  blowing  parameter) 
The  value  of  g  can  be  approximated  using 


g   =   (pu)   St  (2.27) 


From  Reynolds*  analogy  with  unity  Prandtl  number,  the 
Stanton  number  is 


ST   =  -£      =   1 —  (2.28) 

(pu  )p 


where  cf  is  the  local  friction  coefficient.   Thus, 


g   =  —  (2.29) 

UP 


Using  the  Couette  flow  approximation  for  the  boundary-layer 
behavior  with  mass  transfer, 


g   =   g*  £n(l  +  BP)/BP  (.2.30) 


where 


g*   =   lint   (g)  (2.31) 

BP+0 
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In  this  case,  BP  has  been  evaluated  from  the  solution 
of  the  energy  equation  using 

h  -  h. 
BP   =   _P PL.  (2.32) 

bw     fg 

The  wall  shear  stress  has  been  evaluated  using 
Eq.  (2.17)  or  Eq.  (2.20)  and  modified  with  Eq.  (2.30) 

Tbw  =   Tw  £n(1  +  BP)/BP  (2.33) 

where  x   is  the  wall  shear  stress  without  wall  mass  addition, 
w 

The  mass  transfer  conductance  (g)  was  evaluated 
using  Eq.  (2.29).  The  wall  heat  flux  (q")  on  all  solid 
isothermal  boundaries  was  evaluated  using  Reynolds'  analogy 


-qVCh  -  h  )   =  —  (2.34) 

^w7   p    w      u 

p 


Since  the  blowing  rates  were  small  for  the  solid  fuel 
ramjet  (typically  BP  <  2.0),  k  and  e  were  evaluated  using 
Eq.  (.2.4)  and  the  terms  presented  in  Table  II,  which  incor^ 
porate  the  empirical  constants  of  Table  I. 

Five  variables,  namely  u,  v,  k,  e  and  H  Cor)  h  were 
solved  using  Eq.  (.2.4)  in  finite  difference  form.   The  line- 
by-line  iterative  procedure  employed  upwind  differencing 
and  underrelaxation  in  order  to  promote  convergence.   On 
each  radial  line,  the  mass  flow  rate  was  calculated  using  the 
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local  density.   The  "error"  in  mass  flow  (compared  with  the 
summation  of  "mass-in"  at  the  upstream  boundaries)  was  used 
to  uniformly  adjust  the  axial  velocity  over  the  entire  line 
Then,  the  pressure  at  all  downstream  locations  was  adjusted 
to  approximately  correct  for  the  momentum  imbalance  created 
by  the  uniform  axial  velocity  and  a  "pressure  correction" 
equation  was  solved  for  each  cell  on  the  line.   Finally , 
local  cell  velocity  (axial  and  radial)  and  pressure  were 
adjusted  to  satisfy  cell-wise  continuity. 
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III.      ANALYSIS 

A.  INTRODUCTION 

As  previously  discussed,  the  purpose  of  this  study  was 
two-fold;  first,  to  reduce  the  CPU  time  required  to  run  the 
program,  without  considerably  affecting  the  accuracy  of 
results  and  second,  to  make  the  program  output  better  fit  the 
experimental  measurements  without  considerably  increasing  the 
CPU  time.   The  methods  used  in  this  optimization  procedure 
were  the  following: 

a.  Monitor  the  effects  of  the  number  of  traverses  in  the 
tridiagonal  matrix  (TDMA) . 

b.  Monitor  the  effects  of  the  number  of  iterations. 

c.  Monitor  the  effects  of  the  grid  spacing  close  to  the 
solid  fuel  grain  wall,  without  changing  the  number  of 
grid  lines,  in  either  the  x  or  the  r  directions 
(constant  number  of  nodes) . 

d.  Monitor  the  effects  of  changes  in  the  number  of  grid 
lines  in  both  the  x  and  r  directions  (the  number  of 
nodes  was  changed) . 

B.  EFFECT  OF  TDMA  TRAVERSES 

The  number  of  TDMA  traverses  in  the  main  program  is  denoted 
by  the  parameter  NTDMA.   If  there  is  a  predominant  direction 
of  flow,  this  is  arranged  to  be  in  the  x-direction.   Then, 
TDMA  traverses  are  made  on  constant  x  lines,  in  order  to  make 
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the  adjustments  necessary.   Generally,  the  greater  the  number 
of  TDMA  traverses,  the  better  the  accuracy  obtained  but  the 
more  CPU  time  is  consumed. 

C.  EFFECT  OF  NUMBER  OF  ITERATIONS 

The  number  of  iterations  in  the  main  program  is  denoted 
by  the  parameter  LSWEEP.   A  high  number  of  iterations  (or 
sweeps)  is  needed  in  order  to  have  small  residual  values  in 
the  variables  used.   Here  again  the  same  difficulties  as  with 
the  number  of  TDMA  traverses  appeared.   The  larger  the  value 
of  LSWEEP,  the  better  the  accuracy  but  the  greater  the  CPU 
time  required.   It  is  expected  that  for  a  high  number  of 
iterations,  the  increase  in  accuracy  will  be  small  whereas 
the  CPU  time  will  be  excessive. 

D.  EFFECT  OF  GRID  SPACING  CLOSE  TO  THE  SOLID  FUEL  GRAIN  WALL 
In  the  original  program,  the  grid  lines  spacing  close  to 

the  physical  boundaries  was  small.   For  small  grid  spacing 
close  to  the  solid  fuel  grain  wall,  the  non-dimensional  length 
y  at  near  wall  grid  points  was  much  less  than  11.5  and  the 
flow  was  assumed  to  be  laminar.   The  velocity  profile  in  this 
case  is  linear,  given  by 


T 

u   =   — —  6  (.3.1) 

P     ylam 


For  dense  packing  of  grid  lines,  the  flow  may  be  con- 
sidered as  laminar  for  more  than  the  first  node  close  to  the 
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wall.   In  this  case,  the  velocity  profile  will  consist  of 
several  straight  lines  before  turbulence  will  be  encountered, 
as  shown  on  Fig.  (3.1).   This  could  result  in  inaccurate 
predictions.   The  spacing  was  changed  such  that  the  first 
constant-r  grid  line  was  spaced  far  enough  from  the  wall 
in  order  to  have  y   >  11.5.   In  this  case,  the  velocity  pro- 
file will  be  a  logarithmic  curve.   If  this  increased  spacing 
resulted  in  good  predictions,  then  the  number  of  constant-r 
grid  lines  could  be  reduced,  thus  reducing  the  total  CPU 
time. 

E.   EFFECT  OF  CHANGES  IN  THE  NUMBER  OF  GRID  LINES 

Generally,  the  greater  the  number  of  grid  lines  in  both  r 
and  x-directions,  the  better  the  accuracy  expected  but  also 
the  greater  the  storage  as  well  as  the  CPU  time  required  to 
execute  the  program.   Since  it  is  desirable  (but  not  neces- 
sary) to  have  a  maximum  ratio  of  1.5  for  distances  between 
adjacent  grid  lines  in  the  same  direction  and  a  maximum  length 
to  width  ratio  of  10.0  for  any  cell,  it  is  preferable  that 
both  the  number  of  constant-x  and  constant-r  grid  lines  be 
changed  (denoted  by  NX  and  NY  respectively  in  the  main  program) 
In  such  a  case,  if  the  grid  spacing  close  to  the  solid  fuel 
grain  wall  could  be  increased — as  discussed  above — the  number 
of  lines  in  both  directions  could  be  reduced  without  affecting 
the  accuracy  of  predictions.   This  would  result  in  reduced 
requirements  for  storage  and  CPU  time. 
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Figure  3.1 


Velocity  Profile  Close  to  Solid  Fuel  Grain 
Wall 
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IV.       RESULTS 

A.  INTRODUCTION 

All  the  options  discussed  in  the  previous  chapter  have 
been  tried  in  the  modified  CHAMP I0N/2/E/F IX  computer  program. 
The  original  program  had  been  modified  in  the  past  at  the 
Naval  Postgraduate  School  to  account  for  an  aft-mixing 
chamber  [Ref.  2],    and  for  compressibility  and  radiation 
effects  [Ref.  6].   This  program  calculates  the  flow  in  a 
SFRJ,  whose  dimensions  are  shown  on  Fig.  (2.1).   The  input 
values  used  in  the  program  are  shown  in  Appendix  B,  which 
is  a  copy  of  the  BLOCK  DATA  subroutine. 

The  results  of  interest  are  plotted  and  discussed  in 
the  following  sections  of  this  chapter.   The  plots  mainly 
refer  to  two  x-constant  grid  lines  after  the  first  and  second 
sudden  enlargements  at  x  =  202  mm  and  x  =  312  mm,  respectively, 

B.  EFFECT  OF  NUMBER  OF  TRAVERSES 

The  original  program  has  been  modified  to  run  with  three 
different  values  of  TDMA  traverses:   NTDMA  =  3,  NTDMA  =  5 

(originally  used  value)  and  NTDMA  =  7.   Nothing  else  had 
been  changed.   An  attempt  to  run  the  program  with  NTDMA 
greater  than  7  resulted  in  a  repeated  series  of  overflows 
after  a  relatively  high  number  of  iterations  had  taken  place 

(LSWEEP  >  20)  and  the  program  execution  was  terminated. 
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1 .  Temperature  Distribution 

The  temperature  profile  variations  were  negligible 
in  this  case,  as  shown  on  Figs.  (4.1)  and  (4.2).   Thus,  by  using 
a  smaller  number  of  traverses,  the  accuracy  of  temperature 
distribution  prediction  was  not  significantly  affected. 

2.  Effective  Viscosity  Distribution 

Here  again,  the  effect  of  variations  of  effective 
viscosity,  as  shown  on  Figs.  (4.3)  and  (4.4),  were  negligible. 

3.  Axial  Velocity  Distribution 

The  axial  velocity  distribution  after  the  first  and 
second  sudden  enlargements  are  shown  in  Figs.  (4.5)  and 
(4.6),  respectively.   Once  more,  the  number  of  TDMA  traverses 
did  not  affect  the  accuracy. 

C.   EFFECT  OF  NUMBER  OF  ITERATIONS 

The  program  was  run  for  three  different  numbers  of 
iterations;  for  LSWEEP  =  100,  LSWEEP  =  200  (original  valuel 
and  LSWEEP  =  400,  without  any  other  changes. 

1.   y  Distribution 

The  values  of  y  after  the  first  sudden  enlargement 
were  essentially  unchanged.   Some  differences  appeared  after 
the  second  sudden  enlargement  but  these  results  did  not  seem 
to  depend  upon  the  number  of  iterations;  the  peak  value  of 
y+  for  LSWEEP  =  100  was  between  the  peak  values  for  LSWEEP  =  20  0 
and  LSWEEP  =400.   In  any  case,  the  value  of  y  defies 
accurate  prediction  because  of  recirculating  flow;  the  relatively 
high  value  of  y  was  an  indication  of  high  turbulence.   Thus, 
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a  change  in  the  number  of  iterations  did  not  significantly 

affect  the  y  distribution. 
P 

2 .  Temperature  Distribution 

The  differences  in  the  temperature  distributions  after 
the  first  sudden  enlargement  were  negligible,  as  shown  in 
Figure  (4.8).   The  peak  temperature,  2700°K/  which  was  an 
indication  of  flame  location,  moved  almost  to  the  mid-radius 
point.   Then,  a  relatively  steep  reduction  of  temperature 
occurred,  and  at  the  centerline  its  value  was  less  than 
2000°K. 

After  the  second  sudden  enlargement,  the  temperature 
distribution  differences  were  also  negligible,  as  shown  in 
Figure  (4.9) .   The  peak  temperature  now  occurred  at  the  center- 
line  of  the  SFRJ,  which  meant  that  the  enlargement  effectively 
promoted  the  mixing  of  oxygen  and  unburned  fuel . 

3.  Effective  Viscosity  Distribution 

Figure  (4.10)  shows  that  the  variations  in  effective 
viscosity  distributions  after  the  first  sudden  enlargement 
were  almost  negligible.   All  three  curves  coincided  close 
to  the  solid  fuel  grain  wall  and  then  had  a  small  separation 
as  they  approached  the  centerline.   The  same  effect  can  be  seen 
in  Figure  (4.11).   These  small  differences  in  y  ff  close  to 
the  centerline  did  not  really  affect  the  calculations  since, 
as  will  be  seen  later,  the  axial  velocity  distribution  close 
to  the  centerline  was  amost  uniform.   As  a  result,  the  shear 
stress  value  was  also  low  and  was  not  affected  by  the  variations 
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of  y_.cc.   From  the  above,  it  is  seen  that  a  reduction  in 
err 

the  number  of  sweeps  to  one  half  of  the  originally  specified 
value  did  not  significantly  affect  the  effective  viscosity 
distribution. 

4 .   Axial  Velocity  Distribution 

As  for  temperature  and  effective  viscosity  distribu- 
tions/ the  variations  in  the  distribution  of  axial  velocity 
were  also  negligible  for  this  range  of  iterations.   Close  to 
the  solid  fuel  grain  wall,  the  values  of  axial  velocity  were 
the  same  in  all  three  cases.   Close  the  center line  there  was 
a  separation,  with  higher  values  attained  for  higher  numbers 
of  iterations  for  both  of  the  areas  after  the  first  and  the 
second  enlargement  (shown  on  Figures  (4.12)  and  (4.13), 
respectively) . 

Here,  it  was  evident  that  increasing  the  number  of 
iterations  increased  the  accuracy  attained;  although, 
for  this  range  of  sweeps,  the  variation  in  the  predictions 
did  not  exceed  3%. 

D.   EFFECT  OF  GRID  SPACING  CLOSE  TO  THE  WALL 
1 .   y  Distribution 

The  grid  spacing  was  increased  close  to  the  solid 
fuel  grain  wall ,  such  that  turbulent  conditions  would  exist 
at  the  near  wall  grid  points  all  along  the  combustion  length. 
The  values  of  y  have  been  plotted  in  Figure  (4.14)..   Values 
higher  than  11.5  indicate  turbulent  flow,  whereas  lower  values 
indicate  laminar  flow. 
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As  shown  on  Figure  (4.14),  increasing  the  grid 
spacing  close  to  the  wall  resulted  in  the  values  of  y  being 
generally  increased  along  the  region  from  the  first  to  the 
second  sudden  expansion.   After  the  second  sudden  expansion, 
the  values  of  y  that  corresponded  to  increased  grid  spac- 
ing were  generally  lower,  thus  indicating  less  turbulence. 
However,  the  values  were  greater  than  11.5  and  the  flow 
was  handled  as  turbulent  from  the  second  expansion  to  the 
exit. 

2 .   Temperature  Distribution 

The  temperature  distribution  was  dramatically 
affected  by  the  grid  spacing  close  to  the  solid  fuel  grain 
wall  (i.e.,  the  distance  from  the  first  node  to  the  wall) . 
Figure  (4.15)  shows  that  the  original  program  predicted  a 
peak  temperature  rather  far  from  the  wall.   On  the  other 
hand,  the  program  with  increased  spacing  from  the  wall  pre- 
dicted peak  temperatures  close  to  the  wall  (the  larger  the 
grid  spacing,  the  closer  the  distance  of  peak  temperature  to 
the  walll .   A  peak  temperature  farther  from  the  wall  means 
a  thicker  boundary  layer  and,  in  turn,  a  more  fuel  rich 
mixture  ratio  . 

Figure  (4.16)  shows  that  the  original  program  pre- 
dicted complete  mixing  in  the  aft  mixing  chamber .   Peak 
temperature  had  been  predicted  to  be  at  the  center line  of 
the  SFRJ,  whereas  with  increased  spacing  the  peak  temperature 
occurred  just  behind  the  step  wall. 


52 


I      I      t      I      I      I      I      !      I 


I      I      I 


O 

o 
o 
m 


o 

o 

ID 
OJ 


o 
o 
o 

C\J 

IxJ 
CZ 
ID 

\— 

en 

<XI 

• 

o 
o 

LU 

in 


IxJ 


a 
o 
o 


a 
o 
Ln 


o 

OJ 

o 


O 


o 

o 


LP 

o 
o 


a 
o 

a 


a 

<u 

T> 

T3 

CU 

3 

to 

03 

0 

-P 

o 

CO 

u 

cr> 

•H 

c 

PL, 

•H 

o 

0) 

03 

x: 

& 

-p 

CO 

M  T) 

a 

•H 

4-) 

5-1 

M-l 

U 

< 

cu 

c 

x: 

0 

■p 

■H 

■P 

en 

0 

c 

XI 

-H 

•H 

CP 

5-1 

c 

-P 

rti 

w  x; 

•H 

u 

Q 

I* 

cu 

-p 

H 

5-1 

c 

<-| 

3 

cu 

m 

-P 

£  S 

<TJ 

cu 

5-1 

CP 

0) 

cu 

5-1 

X 

a 

m 

-u 

s 

-H 

0) 

c 

0 

E-« 

w 

p 

Ln 


(U 

•H 
Cm 


(W)     SIXU    k 


53 


TD 

- 

•H 

. 

U 

QJ 

o 

cr> 

in 

„, 

LT> 

0 

r-( 

OJ 

CD 

r-l 

-H 

C\J 

V 

C 

U 

to 

•H 

3 

W 

CD 

LU 

to 

c 

cu 

cr 

cy 
u 
u 

-H 

u 

x: 

4J 

O 

o 

o 

c 

a 

0 

OJ 

a: 

M 

to 

-Ul 

cr 

LxJ 
Q_ 

o 

-s~ 

- 

LP 

UJ 

J I L 


!         I 


J I         !         I 


o 

LTJ 

oj 


o 
o 

LO 
OJ 


O 
O 

in 


o 
i_r> 


2 

Q) 
-C 
4J 

c 

a)  o 

en  w 
o 

a  u 
o 

u  tr> 
aj  c: 
en  -h 
u 
<D  rd 
x:  a 

x)  CO 


0) 
4J 


•H 
M 


>4-i  a 

< 

CD 

C  A 

O    4J 
•H 
4-> 

XI 

•H 

U 
4J 
CO 
-H 
Q 

CU 

u 

D 
xJ 
ftJ 
U 
CU 


CT> 

c 

■H 

c 

XJ 

u 


c 

CU 

g 

CU 

tr> 

M 

e  -i 

a  c 


vo 


o 


rvj 
o 


o 


o 
o 

■ 

o 


(wi   Bdidy 


CU 

u 

3 
CP 
-H 
Cm 


54 


3.  Effective  Viscosity  Distribution 

As  shown  on  Figure  (4.17),  the  effective  viscosity- 
distribution  close  to  the  well  was  not  greatly  affected  by 
the  spacing  of  the  grid  lines.   However,  close  to  the  center- 
line  the  effective  viscosity  value  increased  with  increasing 
spacing,  while  the  value  of  yeff  close  to  the  wall  was  slightly 
greater  in  the  original  program. 

The  same  happened  after  the  second  expansion,  as  shown 
on  Figure  (4.18),  although  the  curves  obtained  were  now 
somewhat  different.   Again,  the  effective  viscosity  from  the 
original  program  was  greater  close  to  the  wall,  while,  as 
the  centerline  was  approached,  this  value  dropped  again  below 
the  values  predicted  for  increased  grid  spacing  close  to  the 
wall . 

4 .  Axial  Velocity  Distribution 

As  shown  on  Figures  (4.19)  and  (4.20),  the  original 
program  predicted  much  higher  velocities  than  the  programs 
run  with  increased  grid  spacing  close  to  the  wall .   It  was 
also  noted  that  the  velocity  profile  of  the  original  program 
was  typical  of  laminar  flow,  where  the  ones  for  increased 
spacing  were  typical  of  turbulent  flows  (almost  uniform 
velocity  close  to  the  core) . 

In  this  case  too,  the  same  general  effects  in  the 
distribution  could  be  noted  as  in  those  for  effective 
viscosity.   The  original  program  predicted  higher  values  than 
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the  programs  with  increased  spacing  and  the  .velocity  pro- 
files corresponded  to  laminar  and  turbulent  flows,  respectively 

These  large  variations  in  the  u-velocity  distributions 
shown  on  Figures  (4.19)  and  (4.20)  were  not  surprising  con- 
sidering the  variations  in  the  temperature  distributions. 
Since,  for  increased  grid  spacing,  the  temperature  was  lower 
in  the  core,  the  velocity  was  also  lower  than  that  obtained 
in  the  original  program. 

E.   EFFECT  OF  A  WIDER  GRID  SPACING 

Since  this  case  had  similarities  with  the  previous  one, 
only  one  case  was  run;  with  increased  grid  spacing  in  both 
x  and  r  directions.   The  y  vs.  x  distribution  has  been 
plotted,  as  shown  on  Figure  (4.21) .   It  was  found  that  this 
case  was  not  of  particular  interest.   Along  the  solid  fuel 
grain  wall  y  was  slightly  increased.   This  in  turn  meant 
that  predictions  would  not  be  as  accurate  as  those  obtained 
increasing  the  grid  spacing  close  to  the  wall  only.   So,  no 
further  analysis  took  place  for  this  case. 
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V.   CONCLUSIONS 

From  the  Analysis  of  Chapter  IV  it  can  be  stated  that: 

a.  The  computer  program  was  not  sensitive  to  changes 
in  the  number  of  TDMA  traverses,  NTDMA.   Thus,  a  lower 
value,  such  as  NETDMA  =  3  can  be  used  without  reducing  the 
accuracy  of  predictions. 

b.  As  expected,  the  higher  the  number  of  sweeps,  LSWEEP, 
the  greater  the  accuracy.   However,  if  the  number  of  iterations 
exceeded  200,  the  increase  in  accuracy  was  negligible.   For 
any  number  of  sweeps  between  100  and  20  0,  the  accuracy  of 
prediction  was  found  to  be  very  satisfactory. 

c.  The  program  was  very  sensitive  in  the  grid  spacing 
close  to  the  solid  fuel  grain  wall .   For  increased  grid  spacing 

* 

in  such  a  way  as  to  have  y   >  11.5  for  most  of  the  nodes 
adjacent  to  the  solid  fuel  grain  wall,  the  results  of  the 
program  fit  experimental  data  of  Well  [Ref .  5],  as  shown  on 
Figure  (5.1) .   Also,  this  increased  spacing  close  to  the 
wall  permitted  the' reduction  of  the  number  of  grid  lines 
in  the  r  and,  consequently,  the  x-directions,  with  consider- 
able savings  in  preparation  as  well  as  CPU  time. 
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VI.   RECOMMENDATIONS 

In  view  of  the  results  of  Chapters  IV  and  V,  the  follow- 
ing are  recommended,  regarding  the  execution  of  the  modified 
CHAMP ION/ 2/E/F IX  program. 

a.  The  number  of  TDMA  traverses  should  be  reduced  from 
five  to  three. 

b.  The  number  of  iterations  may  remain  the  same  or  be 
reduced  to  one  half  without  significant  reduction  in  the 
accuracy  of  predictions.   The  use  of  any  number  of  iterations 
greater  than  200  is  not  expected  to  promote  the  accuracy 
significantly . 

c.  While  establishing  the  grid,  the  nodes  which  are 
closer  to  the  solid  fuel  grain  wall  should  be  selected  in 
such  a  way  as  to  result  in  values  greater  than  11.5  (at 
least  for  most  of  them) .   This  will  increase  the  accuracy 

of  predictions  and  will  enable  the  investigator  to  use  fewer 
grid  lines  in  both  the  x  and  y  directions. 
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APPENDIX  A 


C/ASPI#5  JOB  (3238,1085) , 'GENMIX1 /CLASS=C  CHA00 

C/*MAIN  LINES=(12)  CHAOO 

C/  EXEC  FORTXCG,REGION=1024K  CHAOO 

C/FORT.SYSPRIN1  DD  DUMMY  CHAOO 

C/FORT.SYSIN  DD  *  CHAOO 

C  ############################  CHAOO 

C  ANDREA  MH  SVYSHS  AYTO  TO  FAKELLO  TON  XREIAZOMAI  GIA 

TH  THESIS  CHAOO 

C  ################################  CHAOO 

BLOCK  DATA  CHAOO 

COMMON  CHAO  0 

1/CASE51/RIN, IYIN, IYINFl , IXIN, IXINMl ,TAULW (42) , 

XPUSLW (42), CTAULW ,  CHAO  0 

1   CXPLW  CHAOO 
1/CASE52/1XIN2  ,  IYIN2  ,  NXP.l ,  NYRl ,  NXR1M1 ,  NYR1M1 , 

NYRIM2,  CHAOO 

1   CYPTW1 , CTAUT1 , CXPLW2 , CTAUL2 , HWALL2  CHAO  0 
1/CASE1/UINLET , FLOWIN , RP IPE , XP IPE , FXSTEP , HINLET , 

HWALL  CHAOO 
2/DNY/  DYG(42) ,CYV(42) ,FV(42) ,FVNODE(42) ,R(42) , 

RDYG(42) ,RDYV(42)  CHAOO 
2,RSYG(42) ,RSYV(42) ,RV(42) ,RVCB(42) ,RVSQ(42) , 

SYG(42) ,SYGTR(42)        ,  CHAOO 

2,SYV(42)  ,SYVTRV(42)  ,Y  (42)',  YV(42)  CHAOO 
3/DNYONX/AE(4  2)  ,AN(.42)  ,AP(42)  ,AS(42)  ,AW(42)  ,C(42)  , 

D(42) ,DIFE(42)  CHAOO 
3,DIFN(42)  ,DUW(42)  ,DIFW(42) ,DU(42)  ,DV(42) ,EMUE(42)  , 

EMUN(42)  CHAOO 

3,EMUW(42) ,USTAR(42) ,VSTAR(42)  CHAOO 
3,PHIOLD(42) ,RHCE(42)  ,RHON(42)  ,RHOW(42)  ,S  (42)  , 

SPRIME(42)  CHAOO 

3,VQLUME(42) ,CONN(42) ,CONE(42) ,CONW(42)  CHAOO 

4/DNX/DXG(47) ,DXU(47) ,FAREA(.47) ,FU(47)  CHAOO 

4,FUNODE(47) ,KOUNT(.47) ,RDXG(4  7)  CHAOO 
4,RDXU(47) ,RSXG(47) ,RSXU(47) ,STORE(47) ,SXG(47) , 

SXU(47) ,X(47) ,XU(47)  CHAOO 

COMMON  CHAOO 
5/DJPHI/  IEW(IC)  ,ILAST(10)  ,IMON(10)  ,IXNY(10)  , 

IZERC(IO)  CHAOO 
5,KSOLVE(ll)  ,KRS(10)  , RELAX  (10)  ,RSREF(10)  ,RSSUM(10)  , 

TITLE (10)  CHAOO 

COMMON  CHAOO 
6/DO/CCHECK , DP , FLOWPC , FLOWST , FLOWUP , GREAT , IPLRS , 

IPREF,IPRINT  CHAOO 
6 , IWEEEP , IX, IX1NY , IXlNYl , IX2NY2 , IXMON , IXPl , IXPREF , 

IXU,IYMON,IYPREF  CHAOO 


65 


6,JEMU,JH,JLAST,JP,JPP,JRH0,JU,JV,JVP1,KALLED, 

KINPRI,KEHOMU 
6 , KTEST , LABPHI , LSWEEP , NSOLVE , NTRAV , KRAD 
6 , IYF , IYF1 , IYL , IYLK1 , IYLPl , NTDMA , NUMCOL 
6 , N)  , NXMA) , NXM1 , NXM2 , NXYG , NXYP , NXYU , NXYV 
6 , NYNYMA) , NYM1 , NYM2 , P I , RSCHEK , RSMAX , TINY 
6 , IPFRS1 , IPFRS2 , KIN , KEX 
COMMON/PROP , EMUREF , PREFF (10), PRL (10), PRT (10) , 

RHOREF 
COMMON/D2D1/ARSL (47,10), RSLINE (47,10) 
COMMON/D2D2/U(1560) ,V(1560) ,TKEC1600) ,TED(1600) , 
HU600)  ,PP(40)  , 
7P(1444) ,RHO(1600) ,EML(1600) ,RSF(40,40) ,T(1600) 
DIMENSION  F  (14200) 
COMMON 
9/TURB/C1 , 02 , CD , SCRTCD , CD2  5 , ECONST , CTAUTW , CYPTW , 

TAUTW(4  7) ,GENK(42) , 
9YPUSTW(47) ,FACTKE,FACTED,JTKE,JTED,CAPPA,STAN,n^(47) 

HBULKC47) , 
9 RUN , CP , GC , HINR , HWALLR , PREFR , BARM , JT , HOHECK 
9/COMB/AWALL(4  2) ,FLUX(42) ,FLOWALC4  2) ,BP(42) , 

RRC42)  , TEMPI  (.42)  , 
9EMUF , RHCF , HUF , HTR , TREF , STC , HC ,FN2IN ,FOXIN , 
9FMCLFU , FMCLPR , FMCLOX, FMOLN2 
9/CHEMCO/FFU(42) ,FOX(42) ,FN2 (42) , 
9FLFUIT , FLFLOT , FLOXIT , FLOXCT , FLN2 IT , FLN2 ) T , 
9  PCEFU , PCECX, PCEN2 , RRAVE , LINE 
REAL*8  TL(12) 
LOGICAL* 1  LTG(3) 
COMMON 
D/OPLOT/TL 

D,KPLOT(ll)  ,NINT(11)  ,PHI  (40,40)  ,CLC25)  ,XP(47)  , 
SCALEW , SCALEH , LTG 
COMMON/RAD/ TWALL ,TFLAME, SIGMA, EPSIW, IRAD 
E,IRADS,RAD1,FUEL(4  0)  ,QRAD(40)  ,QCON(40)  ,F>!ASF  (40)  , 
VOLL(40) 
DIMENSION  DIFS(42)  ,EMUS  (42)  ,RHOS(.42)  ,CONS  (42) 
EQUIVALENCE  (CONS  (2)  ,CONN(.l)  ) 

EQUIVALENCE  (DIFS (2) ,DIFN (1) ) ,  (EMUS (2) ,EMUN (1) ) 
C      EQUIVALENCE  (RHOS  (2)  ,RHON(-l)  )  ,  (AREAE,AREAW) 
EQUIVALNECE  (RHOS (2) ,RHON(l) ) 
EQUIVALENCE  (F(l) ,U(1) ) 
DIMENSION  A(42),B(42) 
EQUIVALENCE (A CD ,AN(1) ) , (B(l) ,AS (1) ) 

CHAPTER   1 GENERAL, FLOW  PARAMETERS 

DATA  GREAT, TINY, PI/1.E30,   l.E-30,  3.1415926/ 


C 
C 


DATA  RPIPE , XPIPE , UINLET , HINLET , HWALL/ 
10.026924,0.4628,9  2.220,0.0,1.0/ 


DATA  HWALL 2/1.0/ 
DATA  IXIN,IYIN/2,7/ 


CKA00 
CHA0  0 
CHA0  0 
CHA00 
CHA00 
CHA00 

CHA00 
CHA00 

CHAOO 
CHAOO 
CHAOO 
CHAOO 

CHAOO 

CHAOO 
CHAOO 

CHAOO 
CHAOO 
CHAOO 
CHAOO 
CHAOO 
CHAOO 
CHAOO 
CHAOO 
CHAOO 
CHAOO 

CHAOO 
CHAOO 

CHAOO 
CHAOO 
CHAOO 
CHAOO 
CHAOO 
CHAOO 
CHAOO 
CHAOO 
CHAOQ 
CHAOO 
CHAOO 
CHAOO 
CHAOO 
CHAOO 
CHAOO 
CHAOO 
CHAOO 
CHAOO 
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DATA  IXIN,IYIN/2,7/  CHAOO 

C      DATA  KTEST/O/  CHAOO 

DATA  KTEST/1/  CHAOO 

CHAPTER  2 GRID  CHAOO 

DATA  NXMAX,NYMAX/4  0, 40/  CHAOO 

C  CHAOO 

DATA  NX, NY/ 4 0,3 3/  CHAOO 

DATA  KRAD/2/  CHAOO 

C      KRAD=2  MEANS  CYLINDRICAL  POLAR  COORDINATES.  CHAOO 

C      IF  ONE  WISHES  THE  PROGRAM  TO  CALCULATE  THE  GRID,  CHAOO 

C       SET  X(l)  OR  Y(l)  NOT  EQUAL  TO  0.0  CHAOO 

C  CHAOO 

C  CHAOO 

DATA  X/0.0,.01, .023,. 038, .055, .07333, .09166, 

.10999, .12832,  CHAOO 

A. 14665, .1649  8, .18331 ,. 20164, .21997, .2383, .25663, 

.27343, .28468  CHAOO 

B. 29218, .29 718,. 30055, .3028, .304  3, .3053,  CHAOO 
C.3066, .3081, .309 6,. 31185, .31522, .32022 ,. 32772, 

.33892, .35392,  CHAOO 
D. 36892, .38292 ,.39 89 2, .41392 ,. 4289  2, .44192, 

.4543, .4643,  CHAOO 

EO. ,0. ,0. ,0. ,0.,0./  CHAOO 
DATA  Y/0.0, .00061,. 00141, . 00271 ,. 00401, .00531, 

.00631, .00731,  CHAOO 
A. 00861, .01011,. 012235,. 01436, .01586, .01686, 

.01753, .01798, .01828,  CHAOO 

B. 01853, .01873 ,.01888, . 019  00 ,. 01910, .019  25,  CHAOO 

C.  01955, .02015, .02112, .02215, .02315, . 02415 ,. 02515,   CHA01 

D. 0  2565, .026260, .026924,0. ,0.,0.,0.,0.,0.,0.,0.,0./  CHA01 

C  CHA01 

DATA  FXSTEp/1.5/  CHA01 

CHAPTER  3  VARIABLES  CHA01 

C VARIABLES  SOLVED  BY  TDMA  CHA01 

DATA  JU,JV,JTKE,JTED,JH,JPP/1,2,3,4,5,6/  CHA01 

C AUXILIARY  VARIABLES  CHA01 

DATA  JP,JRHO,JEMU,JT,JLAST/7,8,9,10,10/  CHA01 

C —  THE  NUMBER  OF  DEPENDENT  VARIABLES  SOLVED  CHA01 

DATA  NSOLVE/5/  CHA01 

DATA  KSOLVE/9*1,2*0/  CHA01 

DATA  KRS/6*1,4*0/  CHA01 

C  CHA01 

DATA  RELAX/0.3,0.3,0.5,0.5,1. ,1., 0.1. ,01,0.5,1./  CHA01 

DATA  TITLE/1 0*4 HXXXX/  CHA01 

CHAPTER  4 PROPERTY  DATA  CHA01 

DATA  RHOREF,EMUREF/1.403,3.034E-5/  CHA01 

DATA  PRL,PRT/12*1.,   .9,   7*1./  CHA01 

DATA  PREFF/3*1.0,   1.3,   6*1.0/  CHA01 
DATA  C1,C2,CD,CAPPA,ECONST/1.43,1.92,0.09,0.4,9.0/  CHA01 

C      EPSIW  IS  WALL  EMISSIVITY,RAD1  IS  EMP .CONST .FOR  CHA01 

C      EMISSGAS/UNBURNED  FUEL  ON  A  LINE  CHA01 

DATA  SIGMA,EPSIW,TWALL,TFLAME,RAD1/  CHA01 
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15.669E-08, 0.85, 600. , 319  0 .  , 3 .0000/  CHA01 

C      IRAD=0,BP  NOT  MODIFIED  FOR  RADIATION  CHA01 

C      IR/D=1  WITH  RADIATION  CHA01 

DATA  IRAD/0/  CHA01 

C      DATA  IRAD/1/  CHA01 

CHAPTER  5  STARTING  PREPARATIONS  CHA01 

DATA  IXPREF/IYPREF/2/2/  CHA01 

DATA  IXMON,IYMON/22,24/  CHA01 

DATA  KINPRI/0/  CHA01 

DATA  F,DV/14200*0. 0,42*0.0/  CHA01 

DATA  FACTKE,FACTED/0.005,0.0  3/  CHA01 

CHAPTER  6 STEP  CONTROL  CHA01 

CHAPTER  7  BOUNDARY  CONDITIONS  CHA01 

DATA  KIN/KEX/3/l/  CHA01 

CHAPTER  8 ADVANCE  CHA01 

DATA  NTDMA/5/  CHA01 

DATA  RSCHEK/0.001/  CHA01 

CHAPTER  9 COMPLETE  CHA01 

C  CHA01 
C      SUBROUTINE  DENS IT  CALCULATES  VARIABLE  DENSITY 

ON  EACH  LINE.  CHA01 

C      RUN=UNIVERSAL  GAS  CONSTANT  CHA01 

C       CP=C  SUB  P  CHA01 
C       PREFR  IS  THE  REFERENCE  PRES  ADDED  TO  P(I)  TO 

CALCULATE  P  ACTUAL  CHA01 

C  CHA01 

DATA  FLOWAL/42*0.0/  CHA01 

DATA  RUN, OP, GC/ 8 314. 0,1004. 0,1.0/  CHA01 

C  CHA01 

'   DATA  HINP,HWALUR, PREFR, BARM/2. 6570E06, 4. 826E05, 

3.152E05,28.964/  CHA01 
DATA  HTR,FOXIN,FN2IN,STC,HC/2.300E16,0.233, 

0.767, 1.92, 2. 443E07/  CHA01 
DATA  FMOLFU,FMCLPR, FMOLOX,FMOLN2/100 .0,28.97, 

32.0,2  8.0/  CHA01 

DATA  RHOF,EMUF/1199. ,2.43E-0  5/  CHA01 

DATA  TREF/600.0/  CHA01 

C  CHA01 

C      SET  VARIABLE  'LINE'  =  IX  ON  THE  LINE  ON  WHICH 

C       CONSERVATION  CALCULATIONS  DESIRED.  CHA01 

DATA  LINE/29/  CHA01 

C  CHA01 

CHAPTER  10 ADJUST  CHA01 

CHAPTER  11 PRINT  CHA01 

DATA  NUMCOL/10/  CHA01 

DATA  I PLRS, I PRINT/ 2 00, 2 00/  CHA01 

C  ###  CORRECTION  ON  BASIC  PROGRAM  ###  CHA01 

DATA  IPLRS,IPRINT/20 0,200/  CHA01 

CHAPTER  12 DECIDE  CHA01 

DATA  CCHECK/0.0  01/  CHA01 

DATA  LSWEEP.200/  CHA01 
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c######################################## 

C      DATA  LSWEEP/2/ 

C############################################## 

c 

END 
C      PROGRAM  MAIN (OUTPUT, TAPE 6=OUTPUT) 


C 

c 
c 
c 


FOR  CHAMPION  2/E/FIX/CASE50, (HTS7628 ,050 ) 

IDENT  CASE51 
29.7,77 


C/JULY  19  77/ 

C 

C  A  DETAILED  DESCRIPTION  OF  THE  PROGRAM  CAN  BE 

C  FOUND  IN  THE  REPORT   A  GENERAL  COMPUTER 

C  PROGRAM  FOR  TWO-DIMENSIONAL  ELLIPTIC  FLOWS 

C  BY  W.M.  PUN  AND  D.B.  SPALDING,  HTS/76/2, 

C  MECHANICAL  ENGINEERING  DEPARTMENT,  IMPERIAL 

C  COLLEGE,  LONDON 

C 

C 


CHAPTER  0 
C 


0   0   0   0   0 


0   0   DECLARATIONS  0 


COMMON 
1/CASE51/RIN , IYIN , IYINP1 , IXIN , IXINml , TAULW (42), 

XPUSLW(42) ,CTAULW, 
1   CXPLW 
1/CASE52/IXIN2 , IYIN2 ,NXRl ,NYRl ,NXRlMl ,NYRlMl , 

NYR1M2 , 
1   CYPTW1 , CTAUT1 , CXPLW2 , CTAUL2 , HWALL2 
1/CASE1/UINLET , FLOWIN , RPIPE ,XPIPE , FXSTEP , HINLET , 

HWALL 
2/DNY/  DYG(42) ,DYV(42) ,FV(_42) ,FVNQDE(42) ,R(421, 

RDYG(42) ,RDYV(42) 
2,RSYG(42) ,RSYV(42) ,RV(42) ,RVCB(42) ,RVSQ(42) , 

SYGC42)  ,SYGTR(42) 
2,SYV(42)  ,SYVTRV(42)  ,Y(42)  ,YV(42) 
3/DNYONX/AE(42)  ,AN(42)  ,AP(42)  ,AS(.42)„  ,AW(42)  ,C(42).  , 

D(42) ,DIFE(42) 
3,DIFN(42)  ,DUW(42)  ,DIFW(42)  ,DUC42).  ,DV(421  ,EMUE(.421 

EMUNC4  2) 
3,EMUWC42) ,USTAR(42) ,VSTAR(42} 
3,PHIOLD(421  ,RHOEC42)  ,RHON(42)  ,RHOW(_42)  ,SC42).  , 

SPRIME(4  2) 
3, VOLUME  (.42)  ,CONN(42)  ,CONE(42)  ,CONW(42) 
4/DNX/DXGC47) ,DXU(47) ,FAREA(47) ,FU(47) 
4,FUNODE(4  7) ,KOUNT(4  7) ,RDXG(4  7) 
4,RDXU(47)  ,RSXG(47)  ,RSXU(47)  ,STORE(47)  ,SXG(47).  , 

SXU(47) ,X(47) ,XU(47) 


CHA01 

CHA01 

CHA01 

CHA01 

CHA01 

CHA01 

-CCHaOl 

CHA01 

CHA01 

CHA01 

CHA01 

-CCHA01 

CCHA01 

CCHA01 

CCHA01 

CCHA01 

CCHA01 

CCHA01 

CCHA01 

CCHA01 

CCHA01 

CCHA01 

CCHA01 

-CCHA01 

i  CHAO 1 

CHA01 

CHA01 

CHA01 
CHA01 

CHA01 
CHA01 

CHA01 

CHAO  2 

CHAO  2 
CHAO  2 

CHAO  2 

t 

CHAO  2 
CHAO  2 

CHAO  2 
CHAO  2 
CHAO  2 
CHAO  2 

CHAO  2 
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5/DJPHI/  IEW(IO) ,ILAST(10) ,IMON(10) ,IXNY(10) , 

IZERO(IO)  CHA02 
5,KSOLVE(ll) ,KRS(10) ,RELAX(10) ,RSREF(10) ,RSSUM(10) , 

TITLE (10)  CHA0  2 

COMMON  CHA02 
6/DO/CCHECK , DP , FLOWPC , FLOWST , FLOWUP , GREAT , IPLRS , 

IPREF , IPRINT  CHAO  2 
6 , ISWEEP , IX, IX1NY , IXlNYl , IX2NY2 , IXMON , IXPl , IXPREF , 

IXU,IYMON,IYPREF  CHAO 2 
6 , JEMU , JH , JLAST , JP , JPP , JRH 0 , JU , JV , JVPl , KALLED , 

KINPRI , KRHOMU  CHAO  2 
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